function [ s ] = evaluateQuadraturePoint( p, q )
%EVALUATEQUADRATUREPOINT Summary of this function goes here
%   Detailed explanation goes here

r = p - q;

rn = norm(r);

rg = r / rn;

M = [rg(1)*rg(1), rg(1)*rg(2); rg(2)*rg(1), rg(2)*rg(2)];

I = diag(ones(2, 1));

s = -1 / (8 * pi) * (3 * log(rn) * I - M);

end

